Ring structures and mean first passage time in networks 
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In this paper we address the problem of the calculation of the mean first passage time (MFPT) 
on generic graphs. We focus in particular on the mean first passage time on a node s for a random 
walker starting from a generic, unknown, node x. We introduce an approximate scheme of calculation 
which maps the original process in a Markov process in the space of the so-called rings, described 
by a transition matrix of size 0(ln N/ ln(fc) x In iV/ ln(fc)), where N is the size of the graph and 
(k) the average degree in the graph. In this way one has a drastic reduction of degrees of freedom 
with respect to the size N of the transition matrix of the original process, corresponding to an 
extremely- low computational cost. We first apply the method to the Erdos-Renyi random graph for 
which the method allows for almost perfect agreement with numerical simulations. Then we extend 
the approach to the Barabasi- Albert graph, as an example of scale-free graph, for which one obtains 
excellent results. Finally we test the method with two real world graphs, Internet and a network of 
the brain, for which we obtain accurate results. 



PACS numbers: PACS: 



'.75. He, 05.40.Fb, 05.60.Cd 



I. INTRODUCTION 



Modern graph theory starts with the study of Erdos- 
Renyi random networks in 1959 [lj. In more recent times 
it has regained a great amount of attention |2| since it 
has become evident that many different system can be 
described as complex scale free networks i.e assemblies 
of nodes and edges with nontrivial topological proper- 
ties HH. 

In this article we focus on the properties of random 
walks on generic graphs. It is well known that random 
walk is a fundamental process to explore an environ- 
ment [E IE IE IE El j an d recently great attention has 
been devoted to the study of random walk on networks 
(see, for instance, [H [3 HE El EE EE 13 El EE HI ) 
In this process a walker, situated on a given node at time 
t, can be found with probability 1/k on any of the k 
neighbors of that node at time t + 1. 
In particular we are interested in the mean first passage 
time (MFPT) on a node s for a random walker starting 
from a generic, unknown, node x. It is important to note 
here that Noh and Rieger [14| have derived, exploiting 
the properties of Markov chains, an exact formula for 
the MFPT T S j of a random walker between two nodes s 
and j in a generic finite network. 

In this paper, however, we do not trivially average T S j 
over all j : ^ s, a very costly operation, but we use the 
concept of ring (see also |21|). In this perspective we 
study the graph as seen by node s, and partition it in 
rings according to the topological distance of the dif- 
ferent nodes from s (see also |l6t |22|). This allows us 
to map the original Markov problem (of N states) in 
a new Markov chain of drastically reduced dimension 
(0(ln N/ ln(fc) x In NJ ln(fe))) and, as a consequence, to 
calculate MFPT on a generic node s with a reduced com- 



putational cost. On the other hand with the new process 
the identity of the single target node s is lost, and all the 
nodes with the same connectivity (i.e. number of neigh- 
bors) are not distinguishable. 

Our explicit calculation is almost free of approximations 
only for Erdos-Renyi random graphs, for which we obtain 
an excellent agreement between theory and numerical 
simulations. The more disordered scenario of other com- 
plex networks makes the extension of our approach pro- 
gressively more problematic. Nevertheless we find quite 
surprisingly that our approach is able to make very good 
predictions also for other synthetic networks, such as the 
Barabasi- Albert scale free networks |2^| , and at least two 
real world graphs. In all these cases, the considered net- 
works behave, with respect to the property studied, as if 
they were random graph with the same average degree. 
Finally our approach allows us to show that a random 
walker recovers rapidly the degree distribution of the net- 
work it is exploring. 

The paper is organized as follows. In section II the con- 
cept of ring is introduced and the new Markov process 
on which the original problem can be mapped is defined. 
In section III explicit calculations for the case of random 
graphs are performed. It is shown that the description 
of a random graph in term of rings is very accurate, and 
that theoretical predictions for the MFPT are in excel- 
lent agreement with results from simulations. Section IV, 
finally, is devoted to the possible extension of the theory 
to other networks. Notwithstanding the difficulties that 
arise in the analytical extension of the theory, it is shown 
that MFPT of walkers in both artificial and real life net- 
works can be predicted quite accurately with our theory. 
It is also shown that a random walk can be used for the 
reconstruction of the degree distribution of a network. 
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II. THE NEW PROCESS - RINGS 

All the information about a generic graph is contained 
in the adjacency matrix A whose element Ay = 1 if nodes 
i and j are connected, and Ay = otherwise. We shall 
consider here only undirected graphs, i.e. Ay = Aji, 
which do not contain any links connecting a node with 
itself (An — 0, Vi). The degree of a node, k(i), is given by 
k(i) = J2j Ay. Finally, we shall concentrate only on the 
case of connected graphs, i.e. graphs in which each pair 
of nodes i,j, with k(i),k(j) ^ 0, are connected with at 
least one single path. From a random walk point of view 
the matrix A can be interpreted as the N x N symmetric 
transition matrix of the associated Markov process • 

We are interested in the problem of the average MFPT 
on a node s of degree k(s) of a random walker that started 
from a different, unknown, node x. Our idea is mapping 
the original Markov process A on a much smaller process 
B which will be asymmetric and will contain self-loops 
(i.e. Ba ^ 0). More precisely we reduce the N x N 
matrix to a 0(lniV/ln(fc) x lniV/ln(fc}) matrix. 
Given the target node s, we start by subdividin g th e 
entire network in subnetworks, or rings (see also [T6ll22| s ). 
ri , with the following property: 




FIG. 1: Nodes rii per ring: Rings populations are 
shown in this figure. Comparisons of theoretical previsions 
from eq.© and data from simulations for different values of 
k(s) in a single E-R graph of size N = 10 3 and (k) — 6 are 
shown. The fractional error e, defined as the ratio between 
measured and calculated values, is also plotted for each rep- 
resented value of k(s). Theoretical average quantities are in 
excellent agreement with single graph measurements. 



ri = {nodes j | d s j = 1} (1) 

where d S j — dj S is the distance between nodes s and j, 
i.e. the smallest number of links that a random walker 
has to pass to get from j to s. These rings will be the 
states of the new matrix B. Their number, being propor- 
tional to the maximum distance between any two nodes 
in the network , i.e . to the diameter of the network, is 
0(lniV/m(fc)) HHHU (where (k) is the average degree 
of the nodes of the graph). 

Other important quantities are the average number 
m r ^ ri+1 = mj^+i of links that connect all the elements of 
77 with all the elements of r;+i, and the average number 
mi^i of links between nodes belonging to the same rj. We 
have trivially m^-i = toj-i^ and mj.fc = if dik > 1. 
We now have all the elements to define our new process. 
We are no more interested in the exact position of the 
random walker. The relevant information is now the ring 
in which the random walker is. The matrix of this process 
has size (l max + 1) x (l max + 1), where l max is the diam- 
eter of the original graph. The matrix has the following 
structure |33| (for the case l ma x = 6): 
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where fey = TOy/(2fe^7=o mik + 2m ") for * 3-, an d 



bu = 2m ri / (J2k£i=o m ik + 2mii). fey thus represents the 
probability of going from ring i to ring j. By definition of 
rings it is clear that it is not possible to move from a ring 
to a non-adjacent other ring, while it is obviously possible 
to move inside a ring, and in this case the number of links 
must be doubled to take into account that each internal 
link can be passed in two directions. The elements of the 
first row of the matrix are set equal to because we are 
interested in the first passage time in the target node s. 
The probability P$ of going from state i to state j in t 
steps is given by (S*)y. If we set 6qi — F we would allow 
the walker to escape from node s, while feoo = 1 should 
be used if we were interested in the probability that the 
walker reached node s before time t. The probability 
Fk(s)(t) that the first passage on node s occurs at time t 
is then: 

F H*)(t) = E jih {Bt)l ° (3) 
i=i 

where n; is the number of nodes belonging to the ring 
n and each matrix term is weighted with the probability 
that the random walker started in the ring corresponding 
to its row, i.e. ni/(N — 1). 

The average time MFPT r(fc(s)) can be calculated using 
eq.lJSJ as: 

oo 

r(k(s)))=J2 lF Hs)( l ) ( 4 ) 
i=i 
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FIG. 2: Average ring degree: average degrees of nodes be- 
longing to different rings are shown in function of rings' dis- 
tances for an ER random graph (N = 10 4 nodes and (k) = 6) . 
A dependence on ring's distance appears clearly, as predicted 
by theory. Results for different values of k(s) are shown (pre- 
dictions shown only for the case k(s) — 6). 



III. EXPLICIT CALCULATION FOR RANDOM 
GRAPHS 

A. Static 

A random graph is obtained in the following manner: 
given a finite set of isolated N nodes, all the N x (N— 1)/2 
pairs of nodes are considered and a link between two 
nodes is added with probability p. This yields (in the 
limit N — > oo) to Poisson's distribution for the degree k 
of a node: 



P(k) = (5) 

with (k) = p{N — 1). It is clear that such a graph does 
not contain any relevant correlations between nodes and 
degrees, and this will allow us to obtain exact average 
relations for the quantities illustrated in the previous sec- 
tion. 

The first important quantity to calculate is the average 
number n; of elements of 77. It holds 

ni+i = NV-5>*J (l-(l-p)" 1 ) (6) 

where n;+i is calculated as the expected number of nodes 
not belonging to any interior ring that are connected with 
at least a member of r; . Figure Q illustrates that eq. © 
is in excellent agreement with results from simulations. 
Obviously hq — I and n\ = (k). However, if we know 
the degree k(s) of node s we can impose n\ = k(s) and 
calculate the following n;>i in the usual manner. In fact, 
this is the way in which we will use eq.@). 



For the average numbers mij+i of links that connect 
all the elements of r\ with all the elements of and 
ttl\ i of links between nodes belonging to the same rj we 
have: 



mij+i = m yN - ^n k j p (7) 

ni(ni - 1) 
m,i = 2 p 

where the fact that a link between two nodes exists with 
probability p is exploited. As a practical prescription we 
add that when eq. J7J) yields to non-physical mi t i < 
one has to redefine mu = 0. Also these expressions, 
which are crucial for the construction of the B matrix 
of eq. @, give predictions in excellent agreement with 
results from simulations (data not shown). As expected 
we have that, in the limit p — > 0, to/^+i — > while, 
when p 1, rnu+i — * n i x ni+i- Note that from eq. Q 
one has that rii < 1 mij < which has clearly no 
physical meaning. 

Before going o n it is interesting to make a remark. It 
is known [23 . |27| that the nearest neighbors of a generic 
node have particular properties, e.g. an average degree 
different from the (k) of the graph. Our relations are 
able to predict this fact. Combining eq. (JHJ and eq.Q it 
is easy to see that the average degree of the nodes be- 
longing to ri depends on I, being constant (and larger 
than (k)) for low values of I and decreasing rapidly for I 
large enough. Data from simulation reported in Figure |3 
show that this prediction is correct. This agreement is 
not surprising since eq.® and (JJJ are separately in ex- 
cellent agreement with simulations, thus being able to 
predict very accurately the value of the average degree of 
the nodes belonging to the same ring. 

B. Dynamics 

So far, we have shown that predictions of the above 
relations on the static properties of a random graph are 
correct. We now explore the predictions on the diffu- 
sion processes. Figure shows, for different values of 
the degree k(s) of the target node, the comparison of 
the predicted values of i 7 fe( ;s )(<), calculated with eq.© 
and results from simulations. Data from simulations are 
obtained by selecting randomly, for each of the approxi- 
mately N x N runs, one of the nodes of selected degree 
k(s) present in the network to play the role of s, and one 
of the remaining N — 1 nodes to be the starting point 
of the random walker. The agreement between theory 
and simulations is very good. The exponential behav- 
ior, typical of finite ergodic Markov chains, has the form 
f(t) = (l/r)exp(-t/T). 

Figure0]shows that the average MFPT r(k(s)), calcu- 
lated using eq. J3J (with the B matrix built using (eq. © 
and 0), are in good agreement (though slightly smaller) 




FIG. 3: First passage time distributions: FPT prob- 
ability distributions both measured and calculated for an 
ER graph with (k) — 6 are presented for different val- 
ues of k(s). Theoretical predictions are in excellent agree- 
ment with data of random walk on a single graph. Theo- 
retical curves are obtained with eq.Q, and fit the relation 
Fk(a)(t) = (l/T)exp(—t/T) with t obtained with eq.Q. In 
the inset the first part of the distribution is showed in more 
detail. 
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FIG. 4: Mean first passage times: in upper graph MFPT 
both measured and calculated (using eq.lJIJ) are reported for 
an ER graph of size TV = 10 3 with (fc) = 6. Error bars on 
measured values are not visible on the scale of the graph. 
The line r(fe(s)) ~ r(l) x k(s)' 1 is also plotted. It holds 
T(l) s im = 7413 and r(l) ca ; c = 7164 for values obtained re- 
spectively from simulations and from calculation. It can be 
noted that the order of magnitude of r(l) is given by 2M, 
where M is the total number of links in the graph; in our case 
we have 2M ~ (k) x TV = 6000. In lower graph the fractional 
error e, defined as the ratio between simulated and calculated 
MFPTs, is reported. 



with those obtained in simulations. We shall return to 
the origin of the small disagreement between theory and 
simulations at the beginning of the next section. 

The relation, found both in theory and simulations (see 
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FIG. 5: First passage time distribution - Average on 
all target nodes s: The distribution of the FPT obtained 
is shown (ER graph: N — 10 3 , (k) =6). At each run both 
the target node s and the starting node of the random walker 
are randomly chosen. This means that in this distribution 
the degree of the target nodes is no more fixed. The non 
exponential curve is the result of the convolution of several 
exponential curves obtained for fixed k(s). The theoretical 
curve, obtained according to eq.JHJ, is in excellent agreement 
with data from simulation performed on a single graph. The 
fractional error e, defined as the ratio between simulated and 
calculated data, is also plotted. For large times poor statistics 
causes bigger fluctuations of e. 



also T (k( s )) = T (l) x k(s)^ 1 , can be explained with 
elementary qualitative probabilistic arguments. In fact, 
since, as shown in Figure |21 the average degree of the 
nodes in r\ does not depend on k(s) (i.e. on the size n\ 
of ri), also the MFPT on a node of r% is independent 
from m. This means that, while the probability of pass- 
ing from n to s does not depend on n±, a larger n is 
visited more often than a smaller one. Combining these 
observations, it seems plausible guessing that the MFPT 
on a target node s with k(s) > 1 will be 1/k times the 
MFPT of a node s with k = 1, and this behavior is indeed 
observed. 

It is worth noting that both the fc(s) -1 trend and 
the order of magnitude of t(1) can be derived with 
a simple mean field approach 0, ^j. In fact, once 
neglected all the possible correlations in a graph, the 
whole random walk process can be approximated with 
a two state Markov process where the two states cor- 
respond to the walker being at the target node and on 
any other node. Easy calculations shows that the prob- 
ability for the walker to arrive at a node s is given by 
q(s) = q(k(s)) = k(s)/2M, where M is the total number 
of links of the graph. For a fully connected graph this 
relation gives the exact value of r = t(N — 1) = N — 1. 
In a random graph the mean field approach gives bet- 
ter and better results the larger is the mean degree (k). 
For small values of (k), only the order of magnitude of 
t(1) is in fact predicted by this approach. The method 
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based on rings, thought being less simple, is able to make 
more accurate predictions for all values of (k). Just for 
comparison we report here data shown in Figure 0] rela- 
tive to a network of N = 10 3 nodes with (k) = 6 : we 
have r(l) siro /r(l) caic = 1.03 and T (l) sim /2M = 1.24, 
where r(l) s j TO and r(l) ca ; c are MFPT obtained respec- 
tively with simulations and with the ring method. 

All the results discussed above allow us to explain the 
curve presented in Figure [SJ which represents the distri- 
bution P(t) of the MFPT in a graph when both s and 
the starting point of the walker are randomly chosen at 
each run. P(t) can be calculated here as the convolution 
of several exponential FPT distributions F^i s \ (t) corre- 
sponding to the different values of k(s), each weighted 
with the probability of encountering a node of degree 
k(s) in the graph. More precisely, according to eq.JSJ, 
for each time step t every (t) must be weighted with 



Poisson's weights c. 



«fc» = (k) Hs) (k ) 



k(s) 



We have: 



m = E 



„«fe» 



r{j) 



(8) 



This relation can be written in a more compact way ex- 
ploiting the fact that T(k(s)) = r(l) x fc(s)- 1 [gjj. De- 
fined Z = (k) x exp(— t/r(l)), it holds: 
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FIG. 6: First ring fluctuations: Fluctuations related to 
nodes of degree k = 3 are presented for ER (left) and BA 
(right) graphs of size = 1 x 10 J and (k) = 6. On the two 
top figures global fluctuations are analyzed. Histograms of the 
number of nodes vs. the average degree k(r\) of the nodes 
of the first ring r\ for ER graphs (top-left) and BA graphs 
(top-right). Logarithmic scale for the ordinate of BA graph 
must be noted. In the two bottom figures local fluctuations 
are analyzed. Here the histograms represent the number of 
rings in function of the ratio a(ri)/k(n), where a(n) is the 
variance of the degree of the nodes belonging to each r\ for 
ER graphs (top-left) and BA graphs (top-right). It is evident 
the higher degree of fluctuations of the BA graphs. 
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where c is the constant 



IV. EXTENSIONS OF THE THEORY 

In the previous sections we have described a method 
that allows to calculate the average MFPT on a node s 
of a walker that started from a generic other node of the 
graph. We have then obtained exact (average) expres- 
sions for the case of random graphs. Unfortunately, the 
analytical extension of the relations found for this kind 
of graphs to other graphs (such as, for example, scale 
free networks) is difficult. This is due to the fact that 
eq.© and eq.J7|) exploit the knowledge of the rules ac- 
cording to which a random graph is generated. In other 
words the absence of correlations between nodes is the 
main feature those equations are based on. When corre- 
lations are present the calculation of the number of nodes 
of the second ring, , is already very difficult (for finite 
networks) and requires some empirical assumptions |27| . 

In addition there is a more subtle reason that makes 
our method difficult to extend. Given the set of all nodes 
of a graph with a certain degree k, their first rings, al- 
though having the same number of nodes, present two 



kinds of fluctuations. On a global scale, the average de- 
gree k(ri) of the nodes of the first ring does not have 
a unique value, but in general is distributed accord- 
ing to some probability density. On a local scale, on 
the other hand, a single ring is not made by identical 
nodes, and its average degree has a certain variance a. 
In Figure HO we show global and local fluctuations for 
both a random graph and a Barabasi-Albert (BA) scale 
free graph j23(. The preferential attachment rule of the 
Barabasi-Albert network generates a graph with a scale 
free form P(k) ~ k~ c , with c = 3, for the degree distri- 
bution. As it is evident from Figure El BA graphs have 
larger fluctuations than random graphs. 

With our method of rings, described before, the fluc- 
tuations cannot be taken into account. Matrix B J2J is 
in fact defined under both the assumptions of (i) equiva- 
lence of all the nodes with a given k (global homogeneity) 
and (ii) equivalence of all the nodes inside each ring (lo- 
cal homogeneity). The slight disagreement between our 
theory and simulation results present in Figure^are thus 
easily explained in term of local fluctuations of the first 
ring. In fact, as it is easy to demonstrate using Lagrange 
multipliers, the assumed local ordered configuration is 
the most advantageous for a walker that has to reach 
the node s from the first ring. This is then the reason 
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FIG. 7: Mean first passage times vs. the degree of 
the target node for different networks. Continuous 
lines with small filled points (upper graph) are obtained from 
eq.Q, i.e. for random ER graphs, while empty symbols of 
different shapes come from simulations. It is evident the ex- 
cellent agreement between theory and simulations both for ER 
graphs (circles) and, less obviously, for BA graphs (squares). 
Data from real networks are also reported: also in this case 
the agreement with theoretical predictions is good. Error bars 
are not visible on the scale of the graph. In lower graph frac- 
tional errors e, defined as the ratio between measured and 
calculated values, are reported both for real networks and 
BA graph. Dashed lines indicate correspond to e = 1. 



for which our calculated MFPT are always smaller than 
those obtained from simulations. 

Notwithstanding these difficulties in extending our the- 
ory, we found a quite surprisingly result, shown in Fig- 
ure0 given a BA graph with a given average degree (fc), 
the average MFPT for a walker starting from a generic 
node on a node s of degree k(s) is almost equal to the 
corresponding average MFPT of the same random walk 
on a random graph with the same average degree. This 
means that our theory continues to predict very well the 
MFPT (and hence its exponential distribution). It is re- 
markable that the theory predicts well also the MFPT 
on nodes with high degree, which are absent in the cor- 
responding random graph. 

The ability of our theory to predict diffusion processes 
on BA graphs can be due to the modest presence of 
correlations between its nodes. Many properties of real 
networks, in fact, are not reproduced by the BA model. 
One important measure of correlations in a graph is the 
measure of the average degree of the nearest neighbors 
(i.e. of the nodes of the first ring) of vertices of degree 
k, called k nn |2£j- While random and BA graphs have 
a flat k nn , indicating the absence of strong correlations 
among nodes, many real networks exhibit either an as- 
sortative or disassortative behavior. In the first case high 
degree vertices tend to be neighbors to high degree ver- 
tices, while in the second case they have a majority of low 



degree neighbors. Another important measure of corre- 
lation is the clustering coefficient which is proportional 
to the probability that two neighbors of a given node 
are also neighbor of themselves. Again, BA and random 
graphs, in which clustering is very poor, do not reproduce 
the clustering properties of many real networks. 

In order to check how far our theory can predict the 
MFPT on correlated graphs we have performed two sets 
of experiments on real networks. We have considered 
in particular a network of Internet at the level of Au- 
tonomous Systems [29L 130) which exhibits a disassorta- 
tive mixing feature and a recently proposed scale-free 
brain functional network [3ll | which exhibits an assor- 
tative mixing feature as well as a strong clustering coef- 
ficient. 

The results for the MFPT for these two networks, as 
a function of the degree of the target node, are reported 
in Figure Though the agreement between theory and 
simulation is not any more perfect, it remains good. In 
particular we find again the approximate trend r(fc(s)) = 
r(l)xfc(s)- 1 . 

Now we have all the elements to estimate the proba- 
bility of finding a random walker on a node of a given 
degree k. On the one hand, in fact, it seems obvious that 
this probability is related to the fraction f(k) of nodes 
of that degree in the network, while on the other hand 
we now know that the MFPT on such a node is propor- 
tional, on average, to 1/k. It is then reasonable arguing 
that the probability for a random walker being on a node 
of degree k is proportional to kf(k). 
We have tested this hypothesis in an experiment reported 
in Figure |SJ In the experiment a walker has explored a 
BA network and an ER random graph with N = 10 5 
nodes for N time steps. At each time step the degree 
of the visited node was recorded and the normalized his- 
togram of the fraction of time spent on nodes of any de- 
gree is reported in Figure [SJ In the limit of infinite time 
steps this histogram would indicate exactly the probabil- 
ity of finding the walker of a node of a given k. According 
to our hypothesis, this histogram should be described by 
the function P(k)k/{k), where P(k) is the degree distri- 
bution of the considered network, and the Figure shows 
this is in fact the case, already after a relatively small 
number of time steps. 

Finally, it is worth noting that the previous argument can 
be reversed. A walker able to record the degree of each 
node it traverses can be used to determine the degree 
distribution of the network it travels. In fact, if t(k) is 
the fraction of time spent on nodes of degree k it holds 
f(k) oc t(k)/k. The average degree (k) is then trivially 
obtained by requiring the normalization of the estimated 
P(k). 



V. CONCLUSIONS 

In this paper we addressed the problem of the compu- 
tation of the mean first passage time on a selected node 
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FIG. 8: Degree distribution exploration: a random 
walker explores the degree distribution of two networks: a BA 
graph of size N = 10 and (fe) = 6 and an ER random graph 
of size N = 10 J and (fc) = 15. Simulations in which a walker 
travels the networks for N time steps have been performed. 
Empty symbols in figure represent the fraction of time spent 
on a node of degree k. Filled symbols joined by light lines 
are obtained with the relation P(k)k/ (k) , where P(k) is the 
degree distribution of the considered network. Filled points 
fit well the experimental data, and, as it is obvious, a longer 
walking process would allow for better agreement. The prob- 
ability of finding a random walker on a node of a given degree 
k is related to the degree distribution of the graph via the in- 
verse of 1/k MFPT scaling relation (apart from a normalizing 
factor). For BA graphs it holds P(k) oc k~ 3 . The first part 
of the curve P(k)k/(k) presented in Figure is indeed c x fc~ 2 , 
where c is a constant. The region of higher degrees, on the 
other hand, grows linearly due to the fact that in a finite size 
realization the statistics on high degree nodes is poor and the 
degree distribution in this region is flat. We avoided any bin- 
ning in experimental or theoretical data to make clear that 
the walker exploration of the network is not perfect. 



s of random walkers starting from different nodes on a 
generic network. We have introduced a new approximate 
method, based on the concept of rings, which maps the 
original Markov process on another Markov process in a 
much smaller space. This allows for a drastic reduction 
of the computational cost. In the case of ER random 
graphs we have been able to analytically derive all the 
quantities of interest and we have shown that our method 
gives predictions, both for static and dynamic properties, 
in excellent agreement with results found in simulations. 
Even if this new method is promising, analytical results 
arc difficult to obtain for non random graphs. However, 
quite surprisingly, we have found that MFPT calculated 
with our theory for ER graphs are in excellent agreement 
also with simulations of dynamics on BA networks and 
in good agreement with results obtained with random 
walkers on two real networks, thus making our method 
an easy tool to predict MFPT time related quantities in 
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